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Abstract 

We consider the Farey fraction spin chain in an external field h. Using ideas from dynamical 
systems and functional analysis, we show that the free energy / in the vicinity of the second-order 
phase transition is given, exactly, by 

logt 2 t 

Here t = AG'log(2)(l — ^) is a reduced temperature, so that the deviation from the critical point 
is scaled by the Lyapunov exponent of the Gauss map, Xq- It follows that Ag determines the 
amplitude of both the specific heat and susceptibility singularities. To our knowledge, there is 
only one other microscopically defined interacting model for which the free energy near a phase 
transition is known as a function of two variables. 

Our results confirm what was found previously with a cluster approximation, and show that a 
clustering mechanism is in fact responsible for the transition. However, the results disagree in part 
with a renormalisation group treatment. 
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I. INTRODUCTION 



The theory of second-order phase transitions has a long and well-developed history. How- 
ever, for spin models coupled to an external magnetic field there are very few exact micro- 
scopic calculations for the free energy /(/?, h) as a function of both inverse temperature (3 
and magnetic field strength h in the vicinity of such a transition, and (to our knowledge) 
with the exception of the ice-rule models (see Section V. E. of [1]), these are limited to 
mean-field, non-interacting models (e.g. the spherical and Kac-van der Waals models 
or other similarly defined models j^. 

In this work, by using operator techniques, wecalculate the exact free energy /(/3, h) for a 
model with many-body long-range interactions [4]. This model has been investigated previ- 



ously with a cluster approximation [5|, replacing the non-trivial many-body interactions by 
interactions within clusters, and thereby leading to a model very similar to Sj. Our present 
work confirms the validity of that cluster approximation near the critical point. Intriguingly, 
we also find that both critical amplitudes in the free energy scale with a Lyapunov exponent. 

Phase transitions in one-dimensional systems are unusual, essentially because, as long 
as the interactions are of finite range and strength, any putative ordered state at finite 
temperature will be disrupted by thermally induced defects, and a defect in one dimension 
is very effective at destroying order. On the other hand, long range or infinite interactions 
generally make the model ordered at all finite temperatures. Despite this, there are a 
number of examples of one- dimensional systems that exhibit a phase transition. The Farey 
fraction spin chain {4] is one such case, which has attracted interest from both physicists 
and mathematicians (see , llO] and references therein). This model has a phase 

transition at a finite temperature. While the transition is of second-order, it has some 
properties that are usually found with a first-order transition: for external field h = 0, 
the magnetisation jumps from completely saturated, below the transition, to zero above 
it. Despite this unusual behaviour, the model does not violate scaling theory, but rather is 
encompassed limiting case 11 11 1 . 

In some recent work, j^, ll|, this model has been generalised to finite external field, and 
analysed via both renormalisation group methods and with a dynamical system-inspired 
cluster approximation. Neither method is rigorous, and the results are not quite the same. 
Specifically, the dependence of /(/3, h) on h differs by logarithmic factors. This motivates 
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a more rigorous analysis of the model. We find that the cluster picture indeed leads to the 
correct result for the asymptotic form for h), and, in addition, we are able to evaluate 
the constants. Intriguingly, they involve Xq, the Lyapunov exponent of the Gauss map. 
This arises naturally here, since the Gauss map is intimately related to the first-return map 
of the Farey map, which specifies the transfer operator giving the Farey fraction spin chain 
partition function j^. 

In Section [TT] we define the model, first in the standard way using matrices, and then via 
transfer operators. In Section IIIII we derive some operator identities that are necessary for 
our analysis, and in Section [IV] we specify a function space and study the spectral properties 
of the relevant transfer operators. Section |V] describes the connection to the Gauss map. 
Section |Vl] is the heart of our work. Here we use perturbation theory around the critical 
point {l3,h) = (1,0) to find the asymptotic behaviour of the free energy f{j3,h). The key 
to our method is the use of a "cluster operator" , which encodes the behaviour of clusters of 
up and down spins while possessing tractable spectral properties, thus validating the results 
obtained with the cluster approximation of {sj. 



II. THE MODEL 



The Farey fraction spin chain may be constructed, for inverse temperature (3 and magnetic 
field h, via the two matrices 



At 




and A 




(1) 



which correspond to "spin up" and "spin down", respectively. The spin chain partition 
unction comes in various guises, all of which have the same free energy (at least for h = 0, see 



0: 
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12] )■ Here we are considering the generalised Knauf spin chain [6], not the "trace" 



11 



121 for 



model studied in Sj. We make this choice for technical reasons. However, by universality, 
our results are supposed to apply to any of the Farey spin chains (see |6|, 
definitions of the various chains). 
Defining matrix products 

N 

i=l 



(2) 
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where the dependence of Mat on the cxj has been suppressed, and writing a given matrix 
product as = ( " d ) "^^ define the spin chain partition function by 

where x > is a parameter (which does not affect the free energy, as we shall argue below). 
When Mtv begins with A^, c and are neighbouring Farey denominators at level in 



the modified Farey sequence (see ^\ for further details on this connection), whence the 
nomenclature "Farey" for this spin chain model. 
The free energy /(/?, h) is defined via 

- /3/(/3, h) = lim 1 log Z^iP, h- x) . (4) 

N^oo iV 

Alternatively, the partition function can be expressed using transfer operators. In order 
to emphasise the difference between matrices and operators, we denote the latter with script 
letters. We begin by defining the operator 

Cp^H = e-f^'Cl + e^'^jrj, (5) 

where 

= (1^/ (l^) 4/(x) = /(l + x). (6) 



Thus we obtain, as in ^ ^ 

ZN{P,h;x)=C^^M ■ (7) 

This expression indicates that the free energy /(/3, h) is given by the logarithm of the spectral 
radius of on a suitable function space, 

-Pf{P,h;x) = \ogriCp,h) . (8) 

We will return to this point below when we specify the function space used in our analysis. 

There is another notation that is sometimes used in the literature, which we mention for 
completeness and for comparison with previous work on this model. Using a "slash" notation 
which is standard in number theory, the action of a 2 x 2 matrix [M = d) ^ SL2{1j)] on 
a function / is defined via 




(cx + (i)2/^ \cx + d 



Thus ([3]) can be written as 

Zm{(3, h- x) = l(x) I {e-^^A^ + e^^Ayf . (10) 

Note that in order to be consistent with the group structure of S'L2(Z), any addition and 
scalar multiphcation is performed after the matrix action on the function has been computed. 

In the disordered (high-temperature) phase, we expect that there is a leading eigenvalue 
\{P,h) of Cfs^h which satisfies X{P,h) > 1, is non-degenerate, and belongs to the discrete 
spectrum, and that the free energy is given by 

-Pf{P,h;x)=\ogX{P,h) , (11) 

which is independent of x. 

There is an obvious (spin flip) symmetry in our model. Since 



SA^S = A^ with S = I , (12) 




it is natural to define the corresponding operator 5/3 as 

Spfix) = x-'^fil/x) . (13) 

Note that = S and S^^ = S^, i. e. both 5* and are involutions. For the transfer 
operators we find 

Cl = SpCl^Sfi and Ci3,h = SpCfB-hSp . (14) 
The operator Cf^^ has a nice interpretation as the Ruelle-Perron-Frobenius transfer op- 



erator 



13| for the dynamical system given by iteration of the map 



. , , x/(l -x) , < X <1 ^ ^ 

Tix)={ - (15) 

X — 1 , X > 1 , 

and Cfs^h can be viewed as a weighted generalisation. Note that the mapT has the symmetry 
T(l/x) = 1/T(x). This map differs from the Farey map used in, say, [a]. The graph of this 
map is shown in Fig. [H One clearly observes the existence of two marginal neutral fixed 
points, at the origin and infinity. 

Now consider the generating function 

oo 

G{p,h,z;x) = J2^^^N{P,h]x) . (16) 

N=0 







1 



2 

X 
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FIG. 1: The graph of the map y = T{x). 

Examination of G motivates the operator relations discussed below and makes a connection 
with the treatment in [sl- One may rewrite G in terms of the resolvent [1 — zCp^h\^^ of Cp^h 
as 

G{l3,h,z-x) = . (17) 

Equation (fTTll indicates that Zc{l3, h) is equal to the inverse of the spectral radius l/r(£^ /j). 
The free energy is then given as 

(3f{(3,h)=\ogz,{P,h), (18) 

where Zc{(3, h) is the smallest singularity of G{(3, h, z; x) in z on the positive real axis. Thus, 
in principle, we could find the free energy by analysing G. However, it is very difficult to do 
this directly, since Cp^h is not sufficiently well-behaved. In order to make progress we resort 
below to a more nuanced treatment. 
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III. IDENTITIES AND SPECTRAL RELATIONS 



This section introduces a Lemma that is the basis of our analysis. It allows us to avoid 
deahng directly with C^^h, which is difficult to control at the critical point h) — (1, 0). 

To motivate this section, let us consider an arrangement of N spins. We can collect 
adjacent spins of equal orientation into clusters of the form 

n pairs, n > 

>0 >1 >1 >0 

Notice that configurations starting and ending with either spin orientation are included. 
Such an arrangement of spins corresponds uniquely to a particular product of operators 

n pairs, n > 



and £i. 



>0 >1 >1 >0 

and taking a weighted sum over all possible spin configurations of arbitrary length, we find 

[1 - zCp,H]-' = [1 - ze-^'Cl - ze^^'Cl^]-' = 

oo 
n=0 

(19) 

We now introduce the operators 

Ml^^rClil-rCl]-' and A^^,,, = r^fl - r^]-^ (20) 
Notice that as a formal power series in r, 

oo oo 

= (21) 

n=l n=l 

To keep the discussion general, we shall assume in the following that and £^ are bounded 
operators on a Banach space, which will be specified later. 

Then A^jj^^ and -M^^^t- exist whenever the resolvents [1 — r/^j^]"^ and [1 — r£^]^^ exist, i.e. 
for ^ ^i^^p) or T ^ respectively (here cr{A) denotes the spectrum of the operator 

A). This motivates the following identities. 
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Lemma 1. Let z'^ ^ (a{e-f^^Cl) U a{e^^Cj^)y Then 
and 

[1 - zCp,,] = [1 - [1 - -M J,,,,. A^;,,-, J [1 - ^4] . (23) 

Proof. Noting that [l + A^^^^^^] = [l — zCj^]^^ and [1 + -^^ ^e-/3h] = [1 — -^4]"^' calculate 
directly 

[1 + M^^^^^,,] [1 - z/:^.] [1 + A^;,,,-,.] - 1 (24) 

= [1 - zC'^]-'[l - zC,,,][l - zClr' - 1 (25) 

= [1 - ^4]-i(i - zc^,, - [1 - ^4][i - zcimi - zcir' (26) 

= [l-;.4]-M-^4^4][l-^4]-i (27) 

Multiplying ( l22l) by [1 — and [1 — from the left and right, respectively, ( l23|) 
follows. □ 
As above, we have the symmetry 

Ml^ = S.Ml^^^Sp . (29) 

It is helpful to take advantage of this symmetry by defining Mp^r = Aip^Sjs so that 

Note that the rhs is a square when h = 0. Utilising fl30|) and Lemma [H we arrive at the 
following characterisation of eigenvalues and eigenfunctions of Cp^h- 

Proposition 2. Let z'^ i [a{e-^^Cl) U a(e^'^4)) ■ ^/ / 

is an eigenfunction of 

■M.f3^zei3hMi3,ze~i3h With eigenvalue 1, then [1 + ^e-"'']^ ^'^ eigenfunction of C^^h with 
eigenvalue \ = 1/z. Conversely, if g is an eigenfunction of Cp^h with eigenvalue A = 1/z, 
then [1 - zel^^C^p] g is an eigenfunction of M.p^zei^hM.p^ze-i^f^ with eigenvalue 1. 

Proof. If / is an eigenfunction of M-p^zeH'^M-p^ze-fif^ with eigenvalue 1, then by ( !30l) the rhs 
of ( |22l) acting on / is [1 — ■M.'^jj zei^f^'^^ ze-^f^f ~ ^- -^^^ ^'^ assumption on z, the kernels 
of both [1 + A1^_^^^,J = [1 - zef^'^C^p]-^ and [1 + M\^^_fin] = [1 - ^e'^'^^]"^ ^^^O' it 



follows from fl22l) that [1 — zCp^h\g = with g = [1 + -Mp ^^-i3h]f 7^ 0. The second assertion 
follows similarly using ( !23l) . □ 
Proposition [2] motivates the definition of the set 

n^^h = [l/z:zeC\ ({0} U aie-f^Cl) U a(e'^^4)) and 1 G (r{Mp,,eP^Mp,,e-P^)] . 

(31) 

The next proposition relates fljs^h and ai^Cp^h). 
Proposition 3. 

^p,h = \ ({0} U a{e-^^Cl) U a(e^'^4)) . (32) 

Proo/. If z-^ ^ (^(T(e-^'^£j) U ^(e^'^/:^)) , then ((231) implies that [1 - zCp^h]'^ exists if and 
only if [1 — J^^p z^i3hJ^\ ze-i^f^^'"^ exists. But this is equivalent to the definition of Vt^^h, as 
f^/3,hn(a(e-'3'^4)Ua(e^'^4)) =0. 

□ 

Proposition [3] will allow us to study the spectral properties of Cj^^h by analysing the 
spectral properties of ^e/"^•^/3,^e-/5'>• 



IV. THE FUNCTION SPACE 



We now specify the function space on which Cp^h acts and describe some of its spectral 
properties on this space. 

Let Tl = {z E C : 3?2>0} denote the open right half plane and let /7°°(n) denote the 
space of bounded holomorphic functions on 11. Equipped with the norm ||/|| = sup^g^ 1/(^)1 
the space H°°{I\) becomes a Banach space. 

Observe that if is a holomorphic self-map of 11 and w G iir°°(n), then the operator 

: H^ili) ^ H^ili) (33) 

C0/(^) = w{z)f{<P{z)) , (34) 

known as a weighted composition operator (see, for example, [l^), is bounded with operator 
norm ||C^,(^|| = \\w\\. To see this, note that if / G if°°(n), then w ■ f o (p is holomorphic and 
bounded on 11 and 

\\C^J=snp\w{z)f{4>{z))\<\\w\\-\\f\\ (35) 
zen 
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Thus ||C.|i,,<^|| < But ||C^,<^1|| = so = as claimed. 

Before studying the spectral properties of our operators on H°°[Il), we require some more 
notation. We write 

U.J(2) = 1, ,#.l(-) = (l + 2), 

SO that 



(37) 



c 



and Cj^ 



C 



We shall now consider the spectral properties of our operators in more detail. 



(38) 



Proposition 4. 

(i) Cjj is a bounded operator on H°°{Il). Its spectrum is the interval [0,1] with every 
spectral point being an eigenvalue. 

(a) Cj^ is a bounded operator on H°^(n.) provided that 3?/? > 0. If (3 > then Cjj 
for any n G N. 



Proof. For the proof of (i) observe that (f)^ is a holomorphic self-map of 11 and wj^ G H°°[IV) 
so Cj^ is bounded by the argument outlined above. The remaining assertions are proved in 

We now turn to the proof of (ii). Again, since 0^ is a holomorphic self-map of 11 and 
wj^ G /f°°(n) for 3?/? > the operator £^ is bounded by the argument outlined above. 



Suppose now that /3 > 0. For the norm calculation of C^^ observe that 



1 and it follows that 



< 1 for n G N. It remains to show that 



1. Thus 
> 1. In 



order to see this, note that if / G iJ°°(n) is holomorphic at with /(O) = 1, then C}pf is 
also holomorphic at with £^/(0) = 1. Thus £^ 1(0) = 1 for any n G N, and so 
as claimed. 

An immediate consequence of the above is the following. 



> 1 
□ 



I T 

Corollary 5. If f3 > the spectral radii of and Cp are given by 

r(4) = r(4) = l. 



(39) 



The following result will play a crucial role in the study of the spectral properties of Cj3^h- 
Proposition 6. // 3^/3 > then 44 ^-^ compact. 
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Proof. It is not difficult to see that (f)^ o(f)^(Il) = E C : [2; — ||<|}. Thus (f)^ o(f)^ maps 
all of n into a compact subset of 11, and Montel's Theorem jl6l . Chapter 1, Proposition 6] 
implies that Ci ^tot/ii is a compact operator on H°°{Il). But 

and, since o g H°°(Jl), the operator C^T^^i is bounded, so being the product 

of a bounded and a compact operator, is itself compact. □ 
An immediate consequence of the previous proposition is the following estimate for the 
essential spectral radius of Cf^^h (see, for example, [17|, Chap. 1.4]). 

Proposition 7. Let 3fJ/3 > and h E C The essential spectral radius of the operator 
/^p,h = e~^^C^ij + e'^^Cp acting on H°°(Il) is bounded above by e'^*^'^'^^'. 

Proof. Recall that the essential spectral radius of an operator A can be computed as follows 



see, for example, 17|, Chap. I, Thm. 4.10]) 



r^ssiA) = lim ( inf \\A'' - IC\\] ' . (41) 

n— >oo \ /C compact / 

Expanding the n-th power of >C/3,/i, we find 

= e-"'^^4" + e"^^£j" + /C„ (42) 



where /C„ is a sum of 2" — 2 compact operators. In order to see that they are compact, note 



that each them is a product involving the compact operator and bounded operators 



of the form Cr^ and C)^ . Thus 

re..(/:^,/^) < limsup lle-^^'^rj" + e^^'^^^Hi/" (43) 

n— >oo 

and hence, by Proposition |H 

res.(/:Ah) <max{|e-^'^|r(4),|e^''|r(4)} =el^(^'^)l. (44) 

□ 

It turns out that the operator L^^h has a number of interesting spectral properties if 
> and G M, because, in this case, C^^h is a positive operator (see below). Exploiting 

this additional structure requires some more terminology, which we now review (for more 

background see [l8| or [191]). 
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Let H^ijl) = {feH°°{U) : f{z)eRfoTz>0}. This is a real Banach space 
when equipped with the norm inherited from if°°(n) and its canonical complexifi- 



cation equals if°°(n) (see [20|, Lemma 5.2]). In H^{Il) consider the cone K = 
{ / G H^ijV) : f{z) > for z > 0}. Notice that K is closed and reproducing, that is, 
K - K = H^ijl). For f,g e H^iji), we write f < g to mean that g - f e K, and 
this defines a partial order on H^{]1). An operator C on i7]g°(n) is said to be positive (with 
respect to the partial order induced by K) if / > implies Cf > 0, or equivalently, if C 
leaves K invariant. 

I t 

Lemma 8. Let (3 > and G M. Then the operators C^, L^, and Cp^h (^re positive on 
H^{Il) with respect to K. 

Proof. Observe that G K. Thus, f & K implies G K since w^p^z) f {(j)\z)) > for 
z > 0. Hence £^ is positive. By a similar argument, £|j is positive. Finally, C^^h is positive 
since it is a sum of positive operators. □ 
A consequence of positivity is the following lower bound for the spectral radius of Cj^^h- 

Proposition 9. Let > and /i G M. Then the spectral radius of the operator Cp^h = 
e'^^C^l^ + e^'^Cj^ acting on if°°(n) is bounded below by e'^'^'. 

Proof. We start with the case h> 0. Since 

C^ f^l = e-^^w\ + e'^n > , (45) 



the bound r{Cp^h) > follows by [19|, Lemma 9.1], and the positivity of Cp^h- 

Let now /i < 0. For t > define ft{z) := {z + 1)-^^ . Note that ft^K for any t > 0. We 
now claim that for any t > 



Cp,hft > -^^-^ft ■ (46) 
In order to see this note that for z > 

„-f3h e~^^ 1 

where the last inequality follows since {z + t{l + z)) < (1 + t)(l + z) for t, z > 0. Now, as 
before, (146|) and the positivity of Cp^h imply 

r{Cp,,) > (48) 
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by 19|, Lemma 9.1], which in turn yields r[Cp^h) > e by letting t ^ 0. □ 
We now summarise what we know about the spectral properties of Cjs^h- 

Theorem 10. Let /? > and /i G M. For the operator Cp^h acting on H°°{I\) we have the 
hounds 

ress{Cp,H) < el^'^l < r{C(3,h) ■ (49) 

The spectrum of Cp^h in the annulus | z G C : \z\> e'^'*' } coincides with fi^ /j and consists 
of isolated eigenvalues of finite algebraic multiplicity. Moreover, ifr{Cp^h) > ^ess{J^p,h) then 
r{Cp^h) is an eigenvalue of Cp^h. 

Proof. This follows from Proposition [9l Proposition [3] and the definition of the essential 



spectral radius. The last assertion follows from Lemma [8] and 19|, Exercise 8.2]. □ 



We end this section with a number of results on the operators Ai^^^Ai'^^^ which will be 
needed for the perturbative argument below. A short calculation shows that 

oo 

KnKr.f = E n'"^"^'"^ ■ / ° '/'^'"''^^ , (50) 

m,n=l 

with 

^(-'-)(^) = ^ 1 and 0^-'")(^)= '^"^ , • (51) 

^ ^ ^ {nz + mn + 1)2/3 ^ ^ + mn + 1 ^ ^ 

Proposition 11. Let 

D^ = {{P,r^,r^)eC' : 3^/3 > 0, |nl < 1, |r^| < 1 } , (52) 

^2=|(/?,n,rT)GC3 : 3?/?>i |ri|< 1,1x^1 <l| . 
The function {(3, Tj, r^) i-^ A4j^ T^Ai'^p has the following properties: 
(i) on Di it is holomorphic in the operator norm topology; 
(a) on D2 it is continuous in the operator norm topology; 
(Hi) on Di U D2 its values are compact operators. 
Proof. First observe that for 3?/3 > and 2; G n 



(53) 



\nz + mn + 1 



\z 
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so 



hence 



w 



(m,n) 



< 



(mn 



|2K/3 



oo 

< E In 

m,n=l 



m I in 

Inl 



(55) 



(56) 



which means that the series fl5U]) converges in the operator norm topology for any [jS, r^, r^) G 
1^1 U D2. 

Assertions (i) and (ii) now follow by observing that (3 ^ w^^'""^ is holomorphic (and thus 
continuous) for 3fJ/5 > in the norm topology on H°°{Il) for every m, n G N. For the proof 
of (iii) we note that for any fixed 3?/3 > 0, |rj , |r|| < 1 



'Ti-i 



(57) 



so Aijj^^Ai^p^^, is compact by Proposition El The remaining assertion now follows from 
(ii) and the fact that the operator norm limit of compact operators is itself a compact 
operator. □ 

■ ■ 1 IT 

Proposition 12. For /3 > ^7 < r^, < 1, the operator Aij^^^Ai^^^ has a simple leading 
eigenvalue. 

Proof. Fix /3 > |, < r|,r^ < 1. Observe that M^^^^^M^^^ is a transfer operator corre- 
sponding to a real analytic full branch expanding map on [0, 1] with strictly positive weights. 



The proof of 20|, Proposition 4.9] now shows that A^^^^AI^^^ is 1-positive with respect to 
K. Since Ai'^^ ^^Ad^^ is compact by Proposition [TT] the assertion now follows from 18 , 
Theorems 2.5, 2.10 and 2.13] and the fact that the canonical complexification of iJ]^(n) is 



V. DYNAMICAL SYSTEMS 

In this section we provide explicit representations of the operators defined above and 
show the connection to the Gauss map. 

T I 

The operators M.^^., ■M.pr^ ^^"^ ■^P,t have explicit power series expansions in r, given 

by 

n=l ^ ' n=l 

15 



and 



n=l 



'n + x)"^^ \n + x 



(59) 




FIG. 2: The graph of the left-most and right-most branches of the iterated transformations T" 
for n € N. Restricting the iterates as indicated by the solid lines, one obtains the interval map 
y = T{x). 

In order to find an interpretation of these operators as weighted transfer operators asso- 
ciated to interval maps, note that the operators A^jj ^ and TVI^ ^ are given by the collection 
of the left-most and right-most branches, respectively, of the iterated transformations T"' 
for n eN, where T is given in Eqn. (fT5|) . This is indicated in Fig. [2l 

If one restricts the maps as indicated by the solid lines in Fig. [21 one sees that another 
interval map T is formed on R+. Clearly the map T exchanges the intervals (0, 1) and 
(l,oo). 

Recall that the operator Aip^r was defined by composition with as 

Mf3,r = -MJ,.5^ = SpMl, . (60) 

The transformation underlying the operator Ai is obtained from Fig. [2] by composing T 
with S{x) = 1/x. It is easy to see that TS = ST. The graph of TS is shown in Fig. [31 The 
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X 



FIG. 3: The graph of ?/ = TS{x) (soHd hnes), along with the analytic extension of its branches 
(dotted lines). 

associated dynamical system is split into two independent subsystems on the intervals (0, 1) 
and (l,oo), respectively. 

Restricting to act on functions on the unit interval [0, 1] gives precisely the Ruelle- 

Perron-Frobenius transfer operator of the Gauss map x ^ 1/x mod 1, which has branches 

fn{x) = --n. (61) 

X 

Accordingly, Mp^r may be regarded as a generalised transfer operator of the Gauss map, 
with the branches weighted differently. 

This interpretation will be a key ingredient in the perturbative expansion in the next 
section. 

VI. PERTURBATION THEORY 

In this section we come to the final point of our analysis. We employ Proposition [2l 
which connects eigenvalues of Aip^ze/^hJ^f^^^^-ish to eigenvalues of Cp^h to make a perturba- 
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tion expansion around the critical point {P,h) = (1,0). More explicitly, for (3 in a left 
neighbourhood of 1, we will choose z and h so that M-f^^zefhAi/j^ze-f^h has eigenvalue 1. This 
results in an implicit equation for the inverse eigenvalue z{P, h) which by (ITT]) leads to the 
asymptotic form of the free energy h). 

The central object in the perturbative calculation is the operator product 

By Propositions [TT] and [T2l we know that near the critical point, i.e. 2 = 1, /? = 1, and 
h = 0, this is a holomorphic compact operator- valued function that extends continuously to 
this point. Furthermore, at the critical point, "Pi. 1,1 has a simple leading eigenvalue. 
To keep notational overload to a minimum, we write 

f=(n,n) = (ze^^^e-^'^) (63) 

and put Vjs^f = Vp, = M-p^zefi'^Mp^ze-Ph- We further omit any of the variables 

(/3, rj,r|) when they take on their respective value at the critical point, i.e. 2; = 1, /? = 1, 
or /i = 0. For example, we write V = "Pi, 1,1, = 'Pi^zeh,ze-'^ and so on. We use similar 
conventions for other quantities. 
Now we find formally that 

Ch = h with h{x) = - , (64) 

however the function h is not bounded, and therefore not an eigenfunction of the operator 
in the space if°°(n). 

On the other hand, the corresponding equation for V = "Pi, 1,1 is 

Vg = g with g{x) = ^ ] , (65) 

log(2) 1 + X 

and we can check easily that indeed formally h = log(2)(l + Ai\i)g as expected. The 
function g lies in H°°{IV), and is therefore an eigenfunction of V. 

As indicated in the previous section, V is the Perron-Frobenius operator of the second 
iterate of the Gauss map on the unit interval. It follows that the left eigenfunction of V is 
fi = fiL, the Lebesgue measure on [0, 1]. 

Integration of 1/(1 + x) with respect to the Lebesgue measure on [0, 1] gives log(2), which 
motivates the normalisation of g. 

18 



Given the eigenvalue equations 
we shall now solve the equation 

Aftr = 1 (67) 



perturbatively around (3 = 1 and r = (1,1), proceeding as in [22 1. 

By Proposition [11], the compact operator P^g^^ is an analytic function of (3 for 3?/3 > in 
the operator norm topology. Thus, we can expand Vp^i' around /5 = 1 as 

oo 

^/3,r = P.^ + 5^(l-/3)"Pi."^ (68) 

n=l 

Moreover, since by Proposition[T2]the leading eigenvalue .f? is simple, therefore it is analytic 
in (3 and continuous in r. By the same argument, Qis^-f {f^i3,T) is holomorphic in (3 and 
continuous in r with respect to the norm topology on H°°{Il) (the strong dual topology on 
the dual of iJ°°(n)). We thus have expansions analogous to fl68|) for the eigenvalues 
and the left and right eigenfunctions gis^-f^ and /i^^^, respectively. 
We choose the normalisation iip ;fgp ;f = 1. Expanding 

f^TT^/3,Tgi3,T = >^i3,ffJ'Tgi3,f (69) 

and 

fj'P,T'Pp,fgT = >^i3,ff^i3,fgT (70) 

to lowest orders in (1 — f3) and comparing coefficients, we find for the first-order change of 
the eigenvalue 

A?) = f,,V^^^g, , (71) 
which is a standard result of first-order perturbation theory [2^. We therefore have 

Xp^r- = Xr- + (1 - P)firVi''^gr + 0((1 - (3f) (72) 

= A.^ + (1 - P)fiV^''>g[l + 0(1 - + o(l - rj] + 0((1 - (3)') , (73) 

where for the final estimate we have used continuity in r. 

From Vj3 = A^^, where A^/j = 7W/3^i is the transfer operator for the Gauss map, it follows 
that 

^V^^)g = y^M^^^Mg + liMM'^^'^g = 2iiM^^^g , (74) 
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where we have expanded = + (1 - l3)M'-^^ + 0((1 - (3) 



2^ 



By a standard result ISj, this can be expressed in terms of the Lyapunov exponent of 



the associated interval map. Here, one obtains (see e.g. |2J]) the Lyapunov exponent Xq of 
the Gauss map, 



i2 



Therefore Eqn. (|73|1 gives 

A^,.^ = A,^ + 2Ag(1 + 0(1 - + o(l - rj] + 0((1 - pf) . (76) 

Next, we consider the r dependence of X-f = Xi^-f. Using X-ffig-f = fiV^gT^, we rewrite 

fj^Vi'g^ - fiV g^ = {Xf - l)figf . (77) 

Hence 

At = J- H 7 ^ - J- H 7- 7 7 , [ 10 

Kgr) 1 + Kgr - g) 

where we have used the normalisation condition fig = 1. Continuity in f implies that 

A,^ = 1 + ^i{Vr- - V)g[l + 0(1 - r^) + o(l - rj] . (79) 

Combining Eqns. (1671) . (176|) . and (1791) . this implies 

-MPt^-P)^~2Ag(1-/3) (80) 

as /3 — 1 (and, hence, both — 1~ and — > 1^). 

The final step lies in the estimate of — V. We write 

Vf-V = Mr,Mr,-M^ = {Mr,-M)M+M{Mr,-M) + {Mr,-M)iMr,-M) . (81) 

Eqn. (I59l) implies that 

||-Mr--M|| < r/(r) . (82) 

Here 77 (r) = Yl'^=i = Li2(l) — Li2(r), where Li2 denotes the dilogarithm. It follows 
immediately that \\iMr^ - M){Mr-f - M)\\ < ri{Ti)r]{T^). Hence, 

Vf-V = (Mr, - M)M +M{Mr, -M) + 0{r]{Ti)r]{T^)) (83) 

in operator norm. Hence, by (IHOl) . 

- ij{Mr, - M)g + i^iMr, - M)g ~ 2Ag(1 - f3) . (84) 
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An explicit calculation then gives the exact expression 

MA^. - M)g = -^^^ E (85) 

^ ' n=l 

The asymptotic form follows on writing 

" 1 1 

l°g'^ = E^-^-^ + ^(^"')' (86) 

k=l 

where 7 = 0.5772 ... is the Euler-Mascheroni constant. Inserting this in (I85p gives immedi- 
ately 

— ^ 1 — r 1 — r 1 — r 2 1 — r 

n=l 

uniformly for |r| < 1. 

We are now in a position to obtain the asymptotic expansion of the free energy. Inserting 
the leading order asymptotics of fi{M.r — M.)g into equation we arrive at 

[-(1 - n) log(l - n) - (1 - rt) log(l - n)] ~ 2 log(2) Ag(1 - (3) . (88) 

Substituting r| = e^^^^'^'^ and = e^^f~'^'^ and expanding for small / and gives to leading 
order 

2 log(2) Ag(1 -^)-{f + h) log(-(/ + /i)) + {f-h) log(-(/ - /i)) . (89) 

Setting (3c = 1 and C = log(2) Xq = 7r^/6, we see that this is, aside from constants, the 
same as equation (40) in [sl, which was found using a cluster approximation. The analysis 
therein then immediately gives (cf. equation (46) in [sl) 

f-i^.-^T (90) 
logt 2 t 

where the rescaled temperature variable t is given by t = 21og(2)AG(l — (3). Therefore 
the temperature deviation from the critical point is scaled by the Lyapunov exponent of the 
Gauss map, = = 61^(2) ■ ^^^^ that, in addition, (pOl) implies that Xq determines the 
amplitude of both the specific heat C and susceptibility x singularities, with C proportional 
to 1/Xg and x proportional to Ac- 

Also, just as in [s], the asymptotic shape of the phase boundary is given by letting 
— / = \h\ = he in fl89|l . We obtain for the dependence of the critical field strength he to 
leading order 

~ . (91) 

logt 
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VII. DISCUSSION 



It is of interest to discuss the relation of the present treatment to the cluster approxima- 

n 

tion presented in [5|. In that approximation the central quantity is the cluster generating 
function A(/5,r). Here we do not linearise the dynamical map, with the result that the 
function A(/3,r) is replaced by the operator Aiis^r- Conversely, the effect of linearising the 
map on the operator is that it becomes a multiplication operator when acting on constant 
functions. 

Physically, this corresponds to replacing a complicated system with interactions of all 
types (see ^ and references therein) by non-interacting clusters. The significance of our 
work lies in the fact that the behaviour of both models near the critical point is identical, 
thus justifying the cluster approximation. Note that the resulting non-interacting cluster 
model is similar to the ones discussed in j^. 

As mentioned, the renormalisation group result for f{i3,h) found in ll| does not quite 
agree with (|90l) . Specifically, the second term has the form !lJ2e1^ which as t — > 0, is larger 
than the corresponding term in ( l90l) . As discussed in [HI], there does not seem to be any 
consistent way to remove this term in the renormalisation group framework. However, this 
is perhaps not so surprising, since the Farey model is known to have long-range interactions 
(see {4] and references therein), which renders results from a renormalisation group treatment 
questionable. 

Finally, although our results are exact, we indicate which points of our treatment are not 
quite rigorous. The results given in section IIIII are rigorous, but the problem of pointwise 
evaluation is not completely settled, and the spectrum of Vf^^f has not been fully character- 
ized. In particular, the possibility of another leading eigenvalue has not be ruled out. In 
addition, our particular choice of a function space might seem unusual in that it does not 
respect the "spin flip" symmetry of the model. However, it does not seem possible to find a 
function space respecting this symmetry for which the perturbation calculations employed 
are tractable. 
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